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Abstract 

Recent experimental developments showed that the use of the radiation pres- 
sure, induced by a continuous laser wave, to control fluid-fluid interface defor- 
mations at the microscale, represents a very promising alternative to electric 
or magnetic actuation. In this article, we solve numerically the dynamics and 
steady state of the fluid interface under the effects of buoyancy, capillarity, op- 
tical radiation pressure and viscous stress. A precise quantitative validation is 
shown by comparison with experimental data. New results due to the nonlin- 
ear dependence of the optical pressure on the angle of incidence are presented, 
showing different morphologies of the deformed interface going from needle-like 
to finger-like shapes, depending on the refractive index contrast. In the tran- 
sient regime, we show that the viscosity ratio influences the time taken for the 
deformation to reach steady state. 

1 Introduction 

Optical waves are an interesting and promising alternative to pure electric [i] and mag- 
netic fields j3| to produce mechanical stress on fluid interfaces and deform them. 
When focused light propagates through an interface separating two fluids of different 
refractive indices, photon momentum, which linearly depends on the refractive index, ex- 
periences a jump at the interface. This momentum jump induces a radiation pressure on 
the interface that acts towards the fluid of lowest optical index whatever the direction of 
propagation of the optical wave. 

This surprising invariance property was first evidenced in the experiments of Ashkin 
and Dziedzic [3]. Their work was an attempt at giving an experimental answer to 
the Abraham-Minkowski controversy about light momentum's expression in a dielectric. 
The controversy arises due to Abraham's and Minkowski's predictions disagreeing as to 
whether the momentum carried by an electromagnetic field is increased (Minkowski) or 
decreased (Abraham) by the presence of a refractive medium. Ashkin and Dziedzic con- 
cluded that light momentum was consistant with the Minkowski formalism (see [4] for a 
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recent review on light momentum). 

The deflection of the interface towards the less refractive fluid was later confirmed by 
Zhang & Chang [Sj] who used a pulsed laser wave to deform a water droplet and observed, 
immediately after the pulse, an oscillation of the interface leading to an emission of mi- 
crodroplets at the rear of the drop at large pulse energy. More recently, Casner & Delville 
used the soft interface separating two near-critical liquid phases in order to reduce 
by several orders of magnitude the effect of interfacial tension and thus to significantly 
enhance the amplitude of interface deformations. Stationary deformations of large aspect 
ratio were observed [7] as illustrated in Figure [TJ 
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Figure 1: Nonlinear variation of the height h (scaled by the beam radius loq ~ 7.5/im) of the deformed 
interface versus the incident beam power P. The laser wave propagates upwards. The pictures corresponds 
to P = 265, 372 and 478mT/F. Interfacial tension is cr = 5 lO'^iV/m, refractive index contrast is n2 — ni = 
0.012 and density contrast is p^ — p2 — 24.7 kg/m^. 



The radiation pressure effects open interesting perspectives for applications in many 
different fields: 

(i) interface rheology first, with measurements of surface and interfacial tension at small 




scale |8|;|9|, and microrheology characterization |10l: 111 

adaptive optics and holography with light actuation of reconfigurable fluid lenses 
Tsl and interface relief gratings jl3) 

(iii) droplet cavity lasing 

(iv) surface relief micropatterning by extending electrical field approaches jl6| to optics. 
While experiments have highlighted a rich phenomenology, the complete understanding 

of the physics of interface deformation induced by an optical wave is still in its infancy. 
This task is indeed made difficult by the coupling between the shape of the interface 
and the optical radiation pressure distribution as well as the dependence of the interface 
shape on buoyancy, capillarity, viscosity and optical properties of the fluids. While some 
effort has been dedicated to the description of buoyancy effects on equilibrium interface 
deformations of small amplitude [l7| (i.e. linear regime where the deformation's height 
varies linearly with the beam power), very few studies have gone one to describe large 
amplitude deformations (isl; Il9|. 

The aim of this paper is to investigate numerically the statics and dynamics of these large 
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deformations for a wide range of refractive index contrasts, from very soft interfaces (near- 
critical systems cr ~ IQ-'^mN/m) to more usual pairs of fluids (water/air for instance, a ~ 
10 — lOOmN/m ). In essence, it is shown that the underlying physics can be explained by 
the nonlinear dependence of the radiation pressure on Snell-Descartes angles of refraction. 
It is further shown that the transient shape of the deformation exhibits a transition at 
characteristic time which depends on the viscosity ratio. 

Our analysis is organized as follows. The physical model for interface deformation is 
presented in section 2. In section 3, we first validate the computation algorithm through 
a successful comparison with experimental results. Finally an analysis of the refractive 
index contrast and viscous effects is performed showing new important features in both 
the unsteady and steady states in the nonlinear regime of deformation, i.e. when the 
equilibrium hump height has a nonlinear variation with the excitation amplitude. 

2 Governing equations and numerical resolution 
2.1 Physical model 

We investigate the deformation of a initially horizontal stationary liquid-liquid interface, 
induced by the optical radiation pressure due to a continuous Gaussian laser beam of 
power P and radius (also called beam waist) uq. Physical properties of the liquids (de- 
noted 1 for the bottom and 2 for the top) are their refractive indices rii, n2, viscosities /xi, 
/i2 and densities pi, p2- The interfacial tension is denoted by a. The fluids are enclosed 
in a cylindrical cell of radius uq and total height H ^ uq. 

Considering axisymmetry along the beam propagation axis z, cylindrical coordinates 
(Cr, 6^,62) with origin O located at the intersection of the beam axis with the initial 
fiat interface are chosen for this study. A point x is thus referenced by the space coordi- 
nates {r,a,z). The configuration is summarized in figure [21 
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Figure 2: Configuration and notations considered for simulations illustrated using an experimental liquid- 
liquid interface deformation driven by a focused laser wave propagating upward. n2 — n-i — 0.011, 
pi — p2 = 22 kg/w?, a = 3.6 IQ^'' N/m; cjq = 5.3/im ; P — 2AQmW . Si, Sci and Sc2 respectively denote 
the interface and the solid walls in fluids 1 and 2. 
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2.1.1 The optical radiation pressure 

We consider the fluids as two non-absorbing dielectric media with different refractive in- 
dices separated by an interface of arbitrary shape. The indices i and t refer to incidence 
and transmission, and 6i and 6t are, respectively, the Snell incident and transmission an- 
gles (see figure[2]). They can be expressed as 9i = arctan(^) and 9^ = arcsin((l — 6) sin6'j) 
where 6 = is the relative index contrast and z = z{r) is the interface equation. Since 

the photon momentum depends on the refractive index, light momentum is not conserved 
when the beam crosses the interface separating these two dielectrics. The resulting dis- 
continuity in momentum gives birth to the radiation pressure applied to the interface. If 
t and n refer to the tangent and normal directions to the interface at the location where 
the light impinges the interface, z/q represents the optical wave frequency, ho the Planck 
constant and c the light speed in vacuum, the photon momentum balance can be deduced 
from the fact that : 

(i) an incident photon gives the momentum (ni/iot'o/c)(sin6'jt + cos^jii) to the interface, 

(ii) a reflected photon carries the momentum (ni/iot'o/c)(sin^jt — cos6'jn) away from the 
interface and (iii) a transmitted photon picks the momentum {n2hoi'Q/c){sm9tt + cos9tn) 
to the interface. By denoting the number of photons impinging on the interface per 
unit time and unit interfacial area and R{9i, 9t) and T(6'i, 9t) = 1 — R{9i, 9t) the classical 
Fresnel coefficients of reflection and transmission of electromagnetic energy for circularly 
polarized beams are given by j2ol |: 

2(1 — 5) cos6'j cos6't 2(1 — 5) cos cos 

[9,,9t)= ((i_5)cos^^ + coset)2 + {cos9i + {l-5)cos9ty 

We can express the momentum variation dQ, of the incident beam on an interface element 
of area dS during the time dt as rfQ = (iQt + (iQn, i.e.: 

= dQt + (iQn = [rii sin 9i — {Rni sin 9i + Tn2 sin 9t)] '^dSdtt 

+ [ni cos 9i — {—Rni cos 9i + Tn2 cos 9t)]—^dSdtn (2) 

Since ni sin 9i = n2 sin 9t then rfQt = 0. There is no momentum transfer parallel to the 
interface. Consequently, one has 

dQ = rfQ„ = n,cos9Jl + R- ^-^t] ^^^I^dSdtn. (3) 

\ ta,n9t J c 

Classically, the laser intensity / is defined as / = NohoUo, where No is the fiux of photons 
through the beam section. Since the laser wave incidence angle on the interface is 9i, one 
gets N = No cos 9i. We deduce that the optical radiation force acting on the interface per 
unit area dS, when the laser wave propagates from the less to the more refractive fiuid, 
n ^, is normal to the interface and given by 

n-+ r = mcos^^J 1 + i? -^T^^n 4 

\ tan&'i / c 

where 11 ^ is the radiation pressure. 

When the beam is weakly focused, the z dependence of the beam radius can be neglected, 
yielding : 

J(r, z) ^ l[r) = — ^e-'(''/"»)'. (5) 
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Figure 3: Variation of the radiation pressure normalized by its value at normal incidence, versus the 
incidence angle for different values of the relative index contrast S = . 



In figure [HI we have represented 11 ^{6i) normalized by 11 = 0), the radiation pres- 
sure at normal incidence, versus the incidence angle 6i for different values of the relative 
index contrast 5 = ^^^^ = 0.01, 0.02 and 0.4. 

First, we can see that the normalized radiation pressure decreases with 6i and vanishes 
at tangential incidence {6i = 90°). Second, 11 ^{6i) is independent of the incidence angle 
when 6i < 20°, a situation that corresponds to the regime of quasi-normal incidence. For 
6i > 20°, n ^{6i) decreases significantly beyond a certain value of 6i which increases as 
6 decreases. Physically, this illustrates that the radiation pressure, i.e. the momentum 
jump of photons at the interface, is less sensitive to the incidence angle when the two 
fluids tend to refractive index matching (5 — )■ 0). As detailed below, the coupling between 
radiation pressure and interface deformation (tan^j = ^), is the major physical feature 
that explains both the dynamics and statics of the interface deformation. 



2.1.2 Governing equations 

Since the Reynolds number associated with the flows under consideration is always small 
compared to unity {Re = ^^^^ ^ 10^'^, /i ~ 10^^ Pa. s, p ~ lO^fcf^/m^, ooq — 10pm and 
U ~ 10~^m/s a characteristic velocity of the flow), both fluids obey the mass conservation 
and Stokes equations. 

Their dynamics are coupled through stress balance in addition to the continuity of ve- 
locity at the interface denoted Sj. The motion of the interface is described following a 

2a 

Lagrangian approach. Using wq as a reference length, Uq = as a reference velocity 

and piQ = , i = 1, 2 as a reference pressure in fluid i, the axisymmetric dimensionless 
boundary value problem can thus be expressed as 

V.Ui = ; i = 1,2 (6) 
= -Vpi + Aui ; i = l,2 (7) 



5 



:(CTi ■ n - T2 ■ n) = (/s:(r) + n(r) - Bo z)n ; x on Si 



1 + C 

ui = U2 on Si ; ui = on Sci ; U2 = on Sc2 (9) 
— = u(x) on Si (10) 

Ui, Pi are respectively the dimensionless velocity and pressure in each fluid while 
Ti = —pil + (Vui +* Vui) is the dimensionless stress tensor and C = /^i//^2 the viscosity 
Id r— 

ratio. K(r) = -— , is twice the mean curvature of the axisymmetric interface in 

cylindrical coordinates, n(r) = ^11 "^(r) is the dimensionless radiation pressure where 
n~"'' is given by equation (jlj) and Bo = (pi — p2)gojQ/(T is the gravitational Bond number 
which quantifles the relative magnitude of buoyancy to capillary effects for small defor- 
mation amplitude. A detailed investigation dedicated to the gravitational effects, was 
proposed in a previous work, for small deformation amplitudes (itI ; [2iI |. Since Bo is typi- 
cally in the range 10^'^.. 10^^ for fluids commonly encountered, buoyancy can reasonably 
be neglected compared to the other forces and will not be discussed in the present paper. 
In order to quantify effects of the laser wave on the interface deformation, it is convenient 
to deflne the optical to Laplace pressure ratio. This ratio ttq taken at r = (^j = = 0), 
is deflned as: 

TTo = |n(r = 0)1 = — 11) 

vrcwocr {n2 + ni) 

and can be considered as the optical Bond number. 



2.2 Integral formulation and algorithm 

Since solutions to the Stokes problem can be formulated in terms of Green's functions j22|, 
we rewrite the governing equations as a system of integral equations over the boundaries 
of the computational domain. Once boundary conditions on the interface Si and solid 
boundaries Sqi and Sc2 in contact with fluid 1 and 2 are used (see Figure[2]), the two-phase 
Stokes problem can be written in the following compact form: 

1 + C 



u(x) = I V.n{K{ry) + n(rj,) - Bo z{ry))dSy+ (12) 

J Si 

(1-C) f n.'K.MdSy + C I U.(Ti.n)d5y- j U.(T2.n)rf5y. 

J Si J Sci J Sc2 

Here, U and K are Green's kernels for velocity and stress respectively and are given 



by 
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T^/ ,\ 3 ,ddd, ,^ ^, 

where d = x — y, y{ry, Zy) is the integration point. In Eq. ( fT2l) . the flrst term in the right 
hand side describes the flow contribution from interfacial tension, radiation pressure and 
gravity, whereas the second term accounts for shear rates contrast on the interface. This 
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term vanishes when ^ = 1. The third and fourth terms account for shear occurring on 
Sci and Sc2 as a resuh of the no-shp boundary condition. 

Velocities on the interface as well as stress over all the boundaries Sj, Sci and Sc2 
are determined by solving the discrete form of this equation using a Boundary Element 
Method (BEM). Details on the BEM applied to two-phase axisymmetric flow can be found 
in the review by Tanzosh et al. |23|] on the solution of free surface flow using this tech- 
nique. The BEM reveals to be an excellent tool to solve interfacial flow problems with 
high resolution as reported in the analysis of flow involving electric and magnetic fields 
(1| or buoyancy (23; [25 



The solution requires first the discretization of all the boundaries Sj, Sci and Sc2- Due 
to the integral formulation and axial symmetry, the problem is reduced to one dimension 
and only line boundaries need to be discretized. In this work, the mesh is made of con- 
stant boundary elements, i.e. line segments with centered nodes. Azimuthal integration 
of Eq. (fT2|) is performed analytically |26l : [27| reducing Eq. ( IT2|) to a line integration which 
is finally performed using Gauss quadratures j28|. Elliptic integrals resulting from the 
azimuthal integration are evaluated using power series expansions [29.] . 
The fiuid-fiuid interface Sj is parameterized in terms of arc length and is approximated 
by local cubic splines to remesh the interface during its deformation, so that the curva- 
ture can be accurately computed. Distribution and number of points are adapted to the 
shape of the interface, so that the concentration of elements is higher in regions where 
the variation of curvature of the interface is large. The number of mesh points is about 
70 for a typical computation of a small interface deformation. The solid boundaries Sci 
and Sc2 are meshed using about 40 uniformly distributed points. An increase in the mesh 
resolution for the interface and the solid boundaries do not show any change in the results. 
The motion of the interface is followed using the kinematic condition (fTUj) which is dis- 
cretized using an explicit first-order Euler time scheme. A typical computation begins 
with a fiat interface at rest. The laser beam is switched on at t = 0, and the interface 
deforms towards fluid 1 of smallest refractive index. Computation stops when an equilib- 
rium state is reached (v.n|sj- 0). The time step is chosen to be about 20 times smaller 
than the reference time To = 



3 Results 



3.1 Comparison with experiments 

In this section, we compare numerical predictions to experimental data giving both the 
time evolution of the interface and the steady state hump heights and interface shapes for 
various laser illuminations and fluid properties. 

The experimental fluids correspond to a phase-separated near-critical water-in-oil micellar 
mixture contained in a thermo-regulated spectroscopic cell at temperature T, for which 
a very low interfacial tension can be achieved. Details on how the solution is prepared 
can be found in (6|. Above the critical temperature = 308K, two distinct fluid phases 
1 and 2 of different compositions coexist with densities pi > p2 and refractive indices 
Ui < n2- A vertical upward Gaussian beam from a continuous wave Ar+ laser operating at 
wavelength 514.5 nm is then focused at normal incidence onto the flat fluid-fluid interface. 
More exhaustive experimental details on the configuration and protocol as well as on the 
determination of the sam ple properties used here were reported earlier [13; [l9l |. 
As already demonstrated [30|, the thermophysical properties of the fluids can be evaluated 



7 



using asymptotic scaling laws of near critical phenomena |31j . 

Assuming that the coexistence curve is symmetric versus the critical concentration 99c, 
the micellar concentration in each phase ipi (i = 1, 2) can be estimated by: 

^1 = + — , (15) 

Av9 , , 

<^2 = V'c - — , (16) 



with = 0.11 and: 



rj-i rj-i \ 0.325 



= Acp, ( ] , (17) 



with Acpc = 0.5. 

The density of each phase pi, i = 1,2 can be written as a function of ipi: 

Pi = Pmic'Pi + Pconti^ - fi), (18) 

where pmic = 1045 kg.m~^ and Pcont = 850 kg.m~^ are the densities of the micelles and 
continuous phases respectively. 

Because the average distance between two micelles is small compared to the wavelength 
of the laser wave, the mixture can be regarded as homogeneous from the electromagnetic 
point of view. Thus, the relative dielectric permittivity = nf of the mixture is 3^ 



/i2 = ^o( 1-2.5^1, (21) 



/ \ I /I \ V^i(l ^i){^mic ^cont) /-,„n 

ei{(pi) = (fiemic + (1 - f ijf^cont - —/ ^ V- (1^) 

0[^iemic + (1 - '^ijecont) 

along with emic = 1-86 and econt = 2.14, emic and tcont being the relative permittivity of 
the micelles and the continuous oil phase respectively. 

Since concentrations are weak, we use Einstein's relation to estimate the dynamic viscosity 
/ij of each phase: 

Pi = /io (1 + 2.5^] (20) 
2 

with fiQ = 1.269 mPa.s. 

Finally, interfacial tension is written as: 

a = ao l^-^ j , (22) 

with (To = 10"^ N/m. 

Since fluids are assumed to be transparent (the optical absorption is 3 10~^cm~^), the 
action of light on the interface can be considered as a pure mechanical stress, allowing us 
to discard heat dissipation and thermocapillary effects and to consider all liquid properties 
{pi. Pi, rii and cr) as constant and independent of field strength. 

Typical values of the fluids properties are given in table [H 
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T-T,{K) 


a{mN/m) 
1.8 10-4 


n2 - ni 


S 


pi - p2{Kg/m^) 


Ml - P2{Pa.s) 


2 


0.0097 


0.0066 


19 


3.1 10-4 


4 


4.0 10-4 


0.012 


0.0083 


24 


3.9 10-4 


6 


7.0 10-4 


0.014 


0.0095 


27 


4.4 10-4 


8 


1.0 10-3 


0.015 


0.010 


30 


4.9 10-4 


10 


1.3 10-3 


0.016 


0.011 


32 


5.2 10-4 


15 


2.2 10-3 


0.018 


0.013 


36 


6.0 10-4 


20 


3.2 10-3 


0.020 


0.014 


40 


6.5 10-4 



Table 1: Typical values of fluid properties for different values of T — Tc 



3.1.1 Unsteady regime 

One of the first studies of the dynamics of free surface flow induced by the radiation pres- 
sure has been performed by Ostrovskaya et al. j33|. It consisted in solving the unsteady, 
hnearized momentum and mass conservation equations in cyhndrical coordinates with a 
hnearized boundary condition at the free surface. In a recent paper, Wunenburger et al. 
jHI generalized this approach to a liquid-liquid interface where the liquids have indenti- 
cal viscosities. In particular, it was shown that for small deformation amplitudes (where 
h < 1 and ttq 1), the interface dynamics is accurately described by a linear theory 
of overdamped interfacial waves j2]]|. For linearized radiation and Laplace pressures (i.e. 
^ << 1), the following expression of the dimensionless time evolution of the interface 
hump height was predicted to be: 

. ^ 0. ^ ^ f - if (23) 

This prediction has proven to be accurate when compared to experimental data over a 
wide range of fluid properties. 

Note that the equilibrium solution can be easily deduced when t — )■ oo: 

TTn Bo „ , Bo . 

hoo,iin= ^e-E^i—) (24) 

o o 

where Ei(x) is the exponential integral ^-^dk. 

Comparison of numerical predictions with experimental results of the time evolution of 
the hump height for large beam illumination (ttq > 1) is depicted in figure HI Very good 
agreement on hump dynamics and profiles (inset) is illustrated for two sets of experimental 
data (T — Tc = 4.5K and T — T^ = QK) and for different values of the beam illuminations 
ranging from ttq ~ 1.5 to 5.5. 

On the same figure, we have represented the time-dependent solution given by the linear 
model (Eq. [23l) when the optical radiation pressure and capillary terms are linearized. 
This linearization strongly underpredicts both the steady hump height (for T — Tc = 4.5K 
and VTo = 5.5, /ioo,Zm = 4 while h = 7.2) and the time necessary to reach 99% of the steady 
hump height {rgg^m = 40 while rgg = 100). This shows that the complete forms of both 
the radiation pressure and capillary pressure terms are necessary to yield an accurate 
description of the interface dynamics. 
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Figure 4: Time evolution of the interface hump height for ujq — A.Sfim. Experiments (symbols) are 
compared to numerical predictions (continuous lines) and linearized theory hiin{t) (dot-dashed lines). 
(Top)r - = 4.5K, S = ^i^^ = 0.0086 and ( = ^ = 1.37. (Bottom) T -T^ ^ 6K, S ^ 0.0095 and 
C — 1.42. Time t is dimensionless (reference time is tq = ^). The inset shows the interface shape at 
different times. 



3.1.2 Steady state 

In figure El we compare the experimental results for the stationary interface hump height 
h versus ttq to the computed predictions. These data were obtained at T — Tc = lOi^' and 
T — Tc = 3.5K. Experimental profiles of the interface are also compared to the steady 
profiles obtained numerically for ttq ranging from 2 to 4.5. These results clearly show that 
the numerical resolution gives an excellent prediction of the interface hump height and 
profile providing a quatitative validation of our model in both the linear (ttq ~ 2) and 
nonlinear (ttq ~ 4.5) regime. 

Considering the accuracy of our calculations, we now extend the numerical investigation 
to a range of values of the relative index contrast 6 = "S^-j^ii. ^nd viscosity ratio C = — 
that goes far beyond the existing experimental data. This will highlight original features 
for both the steady state and dynamics of the interface. 
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Figure 5: (Top) Numerical and experimental interface steady hump height variation with the optical to 
Laplace pressure ratio ttq. Dashed lines represent the linear prediction hoo,im- (Bottom) Steady state 
interface profiles for different ttq [loq = 5.3/im) (Profiles were arbitrarily shifted in the z-direction). In 
all graphs, symbols and continuous lines represent experimental data and numerical solutions respec- 
tively.(Left) T-Tc = lOK, S = '-^^^ = 0.011. (Right) T-T^ = S.^K, S = 0.008. 



3.2 Optical effects 

The near-critical fluid phases used in the experiments had a relative index contrast varying 
from 6 = 0.0066 for T — Tc = 2K to 6 = 0.015 for T — Tc = 25K featuring a very narrow 
range of variation of this parameter. Experimentally, it is possible to reach steady defor- 
mations with comparable aspect ratios using other couples of fluids such as Octane/ Water 
with Span 80 or SDS as surfactants (a ~ 10~^mN/m and 6 ~ 0.06). Higher relative in- 
dex contrasts 6 ~ 0.35, (Water+Lung surfactants) /Air interface (a ~ lOmN/m) can be 
achieved, altough a beam power of approximately lOW becomes necessary. 
Given these experimental constraints, we investigate the role of 6 in the range 10~^ < 6 < 
0.4 by means of numerical simulations in order to understand its effects on the deforma- 
tion in the nonlinear regime. 

In flgureini we have represented the variation of the steady state interface hump height, 
h as a function of ttq for several values of 6 ranging from 0.01 to 0.4. 
As shown above, for sufficiently small values of the optical illumination (ttq < 2), h varies 
linearly with ttq and has no dependence on 6. 

In the nonlinear regime (ttq > 2), h strongly depends on 6. For any given value of ttq, h 
is larger when 6 decreases. This is due to the fact that, for given ttq and 6i, n(6'j) is a 
decreasing function of 6 (see figure E]) • 

The inset of figure [H] illustrates this statement. It shows that the value of ttq needed to 
reach /i = 6 as a function of 6 decreases when decreasing 6. However, we can also notice 
a saturation of ttq = f{S). Indeed, when 6 is getting very small, the radiation pressure 
becomes quasi insensitive to the incidence angles 9i (i.e. to the local slope ^) except 
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Figure 6: Interface hump height variation with ttq for different relative index contrasts S — ■ The 

square symbols are experimental data corresponding to (5 = 0.01. (Inset) Variation of tto needed to reach 
ft. = 6, as a function of b. 

when Qi is very close to 90°. 

In figure [7] (a-b-c) we have represented the interface steady profile, for three values of 
5: 10~^, 0.01 and 0.4. The corresponding values of ttq were chosen in order to obtain the 
same interface hump height /i ~ 10. It is noteworthy that the shape of the deformation in 
the nonlinear regime strongly depends on b. At relatively large values of 5, (e.g. 5 = 0.4 
in figure [7^), the center part of the interface adopts a needle-like shape with a sharp tip. 
Decreasing 5 (e.g. 6 = 0.1, see figure [7b) widens and rounds off the tip, showing a bell 
shape, while the base of the deformation shrinks. More surprisingly, when the two fluids 
have almost the same refractive properties, a stable finger, nearly cylindrical with a round 
tip, is obtained (figure [7b). This shape evolution corresponds to what is observed on the 
two pictures superimposed on figure [7| and extracted from experiments performed with 
near-critical phases. Figure [TJl, obtained at T — Tc = 20K {6 ~ 0.014) and ttq ~ 5.5 
shows that the interface adopts a bell shape while figure [7^, obtained at T — = 2K 
[6 ^ 0.0066) and ttq ~ 9, shows a nearly cylindrical finger. 

These observations can be explained qualitatively by considering the balance of radiation 
and Laplace pressures. The finger shape {6 — )■ 0) corresponds to a case where the ra- 
diation pressure is almost independent of the incident angle (see figure [3]). This means 
that n(r) TToC"^'' . Therefore, as only the capillary force balances the optical radiation 
pressure at steady state (we neglect the weak infiuence of gravity), we have /t(r) ~ n(r). 
This is confirmed on figure [3-h where the curvature profile shows an excellent agreement 
with the gaussian profile (see figure [7h). 

The other shapes (5 = 0.1 and S = 0.4) show that when the radiation pressure strongly 
depends on the incidence angle, the profile tends to a needle shape as the maximum of 
n(r) becomes more and more localized toward the tip (r = and 6i = 0). 
Consequently, one can tune the morphology of stationary nonlinear deformations by sim- 
ply changing the refractive index contrast. 
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Figure 7: Bottom (a-b-c): Interface steady profiles for S — 0.4, 0.01 and 10 ^. ttq is chosen to keep 
h ~ 10. Bo — 0.003. Central pictures (d-e) were obtained from experiments with Bo — 0.003. Left 
picture (d): T-Tc = 20K, S ~ 0.014, ttq ~ 5.5. Right picture (e): T-Tc = 2K, 6 ~ 0.0066 and ttq - 9. 
Top (f-g-h): Curvature of steady profile associated to the deformations. The line with (+) symbols is a 
gaussian profile. 

In the following, we show that the dynamics of the deformation has also a strong depen- 
dence on the relative index contrast 6. 

In figure El the time evolution of the hump height h{t) is plotted for 5 ranging from 
10~^ to 0.4. We note that the transient regime is significantly longer when 6 is closer to 
zero, i.e. when the fluids tend to match optically. This is explained in the left inset of this 
figure that shows the time variation of the radiation pressure integrated over the interface 

R R 

I Il{r,t)ndS\ = 2tt Jq"° Il{r,t)rdr reduced by its initial value 271 Jq"° n(r, t = 0~^)rdr. 
This ratio, called F, simply represents the time variation of the net optical force applied 
to the interface during the transient deformation. We can notice that this force quickly 
reaches its steady state value when S = 0.4 whereas the transient is much longer for 
6 = 0.01. This observation can be explained qualitatively by considering the decrease 
of the radiation pressure with the incidence angle 6i. As the deformation increases, 6i 
increases along the interface, therefore the associated radiation pressure applying on the 
growing hump decreases locally for S = 0.4 while it remains quasi unchanged for 6 = 0.01. 
The characteristic relaxation time of the optical force is smaller when increasing 6. As this 
force is the source of momentum transfered to the interface, the dynamics of the interface 
follows a similar behaviour with 6. 

The right inset shows a characteristic time Tqq chosen arbitrarely at 60% of the equilibrium 
hump height as a function of 6. We can clearly see that Tqq increases when 6 decreases, 
before reaching a saturation. This saturation is related to the weak dependence of n(r) 
on 6i at small 6. 
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Figure 8: Time evolution of the interface hump height for different values of i5. ttq is chosen to keep /i ~ 6. 
Bo = 0.005. The left inset shows the time variation of the net radiation force acting on the interface F{t) 
for S = 0.4 and 6 = 0.01 (C = ^ = !)■ The right inset shows the characteristic time tqq taken for the 
hump height to reach 60% of its final steady value plotted against 6. 



3.3 Viscous effects 

By analogy with tlie analysis of the effects of refraction indices on the deformation, we 
now investigate the effect of the viscous ratio C = — on the transient behaviour of the 
deformation towards equilibrium. Equation ([H]) shows that the transient deformation de- 
pends on viscosities only through their ratio (. 

In the case of the near-critical phases used in the experiments reported earlier, ( was 
varying from 1.27 for T — Tc = 2K to = 1.76 for T — Tc = 25K. As ( was close to unity, 
no particular behaviour was observed in the transient regime. 

It is however possible to imagine experiments with larger values of (. Using viscous 
Glucose- Water/Hexadecane or Glucose- Water/Toluene and appropriate surfactants would 
enable to obtain a viscosity ratio of ~ 150 with 6 ~ 0.06 and cr ~ IQ-^mN/m still al- 
lowing large scale deformations using a continuous laser beam. A theoretical study would 
thus be a very useful predictive tool for the design of realistic future experiments. For 
this reason, we investigate the effect of the viscosity ratio in the two limit cases ( = 200 
and C = 1/200. 

Figure [9] shows the evolution of the interface profile at different times for these two 
values of (, keeping 6 = 0.1, Bo = 0.01 and ttq = 15. 

In the case C. = we can observe a sharp transition from a bell shape (t = 0.65), 
characteristic of a small deformation amplitude, to a needle shape (t = 4.7). For ( = 200, 
we observe that the transition towards the needle shape occurs later, the tip of the in- 
terface showing a rounded shape between t = 4.7 and t = 7.9. This shape transition can 
be further characterized by analyzing the time dependence of the reduced curvature of 
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Figure 9: (Top) Interface profiles at different times for two values of the viscosity ratio C = ^ = 2M 
(left) and ( — 200 (right). S = 0.1 and Bo = 0.01. Profiles were shifted vertically for clarity. (Bottom) 



Time evolutions of the corresponding reduced curvature of the curvature — = — — ^— ^. 



^(2)(t) ^^^^^ 
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the curvature at the tip, — = ^ = — ^ — - — shown in figure [HI 

K{t) K{r = 0,t) j-§\ 



> (t) r dr y dr 

(Pz_\ • K 

dr2 l''-=0 

being an indicator of the spatial variation of the curvature. In fact, a needle shape has a 
sharper curvature variation with r than a rounded shape (see figure [Tff) and (h)), which 
means that is larger (in magnitude) for the needle shape. We observe in figure 
that for a given value of '^^K(t)^ initially increases until reaching a maximum at time 
(ti/200 = 3.3, T200 = 9.5) before decreasing toward the steady state value . This maximum 
is a signature of the change of the shape evolution that switches from a round tip to a 
sharp one corresponding to a needle shape. 



4 Conclusion 

While quantitatively validated with experimental data in the nonlinear regime of deforma- 
tion, our numerical resolution, based on a Boundary Element Method, showed that new 
morphologies of a soft fiuid-fiuid interface deformed by a continuous laser wave emerge 
when fiuids become contrasted. In the nonlinear regime, where the radiation pressure 
and the height of the interface strongly depends on the relative refractive index contrast 
6, we showed that the interface shape turns from a needle to a nearly-cylindrical finger 
shape. In the transient regime, we predicted that the characteristic time of the defor- 
mation increases when decreasing 6 before reaching a saturation. These results show the 
strong nonlinear coupling between the radiation pressure and the interface deformation 
in both transient and steady state. It was shown that the physical feature explaining 
these results lies in the dependence of the radiation pressure on the incidence angle. An 
original experimental evidence of this dependence can be achieved when the beam axis is 
not perpendicular to the initial fiat interface. In this configuration, the interface hump is 
attracted in the direction of light propagation in the nonlinear regime of deformation (see 
figure [lO]). 

Finally, the infiuence of the viscosity ratio ( has been investigated for large scale defor- 
mations. We showed that the transition time of the interface from its initial bell shape to 
the steady shape increases with (. 

All these results show that the morphology of nonlinear interface deformations driven by 
the radiation pressure of a laser wave is even more rich than expected, opening the route 
towards future original experiments. Further developments on nonlinear effects, where a 
feedback coupling between deformation and propagation emerges, will extend deformation 
shapes to self-adapted ones as for liquid optical fibres [34] and nipple-like shapes observed 
experimentally |35|. They both involve light guiding within the deformation that modifies 
in turn the radiation pressure along the structure. This thorough analysis of optohydro- 
dynamics provides the general frame to predict and anticipate further developments of 
contactless interface manipulation at the micrometer scale. 
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Figure 10: Interface deformation with a tilted laser beam propagating upwards. Top : P = IW . Bottom 
P = 1.7W. n\ — n2 = 0.016 and a = 1.4 10~^N/m. When the beam is tilted, incidence angles on one side 
of the deformation (here the left one) are larger than in the other side. When they are sufficiently small 
to consider radiation pressure as constant all over the interface, the deformation remains symmetric (Top 
picture). When the beam power is increased, the growth of the deformation makes radiation pressure 
sensitive to incidence angles and then more efficient on the right than on the left. The deformation is 
thus attracted in the beam propagation direction and becomes non-axisymmetric (Bottom picture). 
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